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ABSTRACT 


The  maximum-flow  network-interdiction  problem  (MXFI)  arises  when  an 
interdictor,  using  limited  interdiction  resources,  wishes  to  restrict  an  adversary’s  use  of  a 
capacitated  network. 

MXFI  is  not  easy  to  solve  when  converted  to  a  binary  integer  program.  Derbes 
(1997)  uses  Lagrangian  relaxation  to  solve  the  problem,  at  least  approximately,  for  a  single 
value  of  available  resource,  R.  Bingol  (2001)  extends  this  technique  to  solve  MXFI 
approximately  for  all  integer  values  of  R  in  a  specified  range.  But,  “problematic  /?- values” 
with  substantial  optimality  gaps  do  arise.  We  reduce  optimality  gaps  in  two  ways.  First, 
we  find  the  best  Lagrangian  multiplier  for  problematic  /?- values  by  following  the  slope  of 
the  Lagrangian  function.  Secondly,  we  apply  a  limited-enumeration  branch-and-bound 
algorithm. 

We  test  our  algorithms  on  six  different  test  networks  with  up  to  402  nodes  and  1826 
arcs.  The  algorithms  are  coded  in  Java  1.3  and  run  on  a  533  MHz  Pentium  III  computer. 
The  first  technique  takes  at  most  39.3  seconds  for  any  problem,  and  for  one  instance,  it 
solves  8  of  the  problem’s  15  problematic  i?-values.  For  that  problem,  the  second  technique 
solves  four  of  the  remaining  problematic  /?- values,  but  run  time  increases  by  two  orders  of 
magnitude. 
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EXECUTIVE  SUMMARY 


This  thesis  is  concerned  with  solving  or  approximately  solving  a  maximum-flow 
network-interdiction  problem  (MXFI).  This  problem  arises  when  an  interdictor,  using 
limited  interdiction  resources,  wishes  to  restrict  an  adversary’s  use  of  a  capacitated 
network.  In  particular,  the  adversary,  or  “network  user”  seeks  to  maximize  flow  through 
the  network,  and  the  interdictor  seeks  to  minimize  that  maximum  flow  by  destroying  arcs. 

MXFI  is  not  easy  to  solve  when  converted  to  a  binary  integer  program.  Derbes 
(1997)  uses  Lagrangian  relaxation  to  attempt  to  solve  this  problem  for  a  single  value  of 
interdiction  resource,  R.  Bingol  (2001)  extends  this  technique  to  solve  MXFI  for  all  integer 
values  of  R  in  some  specified  range  assuming  that  exactly  one  unit  of  resource  is  required 
to  interdict  an  arc.  However,  about  25%  of  these  values  are  “problematic”  in  tests,  i.e., 
positive  optimality  gaps  arise. 

In  our  study,  we  first  focus  on  decreasing  the  optimality  gaps  by  finding  the  best 
Lagrangian  multiplier.  We  simply  follow  the  slope  of  the  Lagrangian  function  to 
accomplish  this.  For  most  test  problems,  this  method  decreases  the  gap  for  problematic  R- 
values  by  approximately  50%.  This  method  helps  find  better  feasible  solutions  and 
sometimes  prove  optimality.  But,  some  large  optimality  gaps  may  remain. 

We  tackle  the  remaining  problematic  /Lvalues  with  a  restricted  branch- and-bound 
technique,  which  has  only  one  level  of  enumeration  below  the  root.  This  technique  creates 
a  strict  partition  of  the  feasible  region  and  the  minimum  lower  bound  from  among  these 
partitions  is  a  lower  bound  on  the  optimal  solution  value.  In  our  tests,  we  solve 
approximately  70%  of  the  problematic  /Lvalues  that  remain  after  applying  the  first 
technique.  However,  branch-and-bound  increases  run  times  dramatically. 

Both  procedures  are  written  and  compiled  using  the  Java  1.3  programming 
language.  All  tests  are  performed  on  a  personal  computer  with  a  533  MHz  Pentium  III 
processor  and  192  MB  of  RAM,  running  under  the  Windows  XP  operating  system.  For  test 
purposes,  we  have  used  6  different  grid  networks.  The  smallest  one  has  27  nodes  and  86 
arcs  and  the  largest  one  has  402  nodes  and  1826  arcs. 
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Applying  both  techniques  together,  we  have  been  able  to  either  solve  problematic 
/?- values  exactly  or  improve  their  optimality  gaps  for  most  cases.  In  fact,  85%  of  the 
problematic  /?- values  are  solved  optimally  in  our  tests.  One  problem  with  a  large 
optimality  gap  does  not  improve,  however. 


I.  INTRODUCTION 


This  thesis  develops  a  method  based  on  Lagrangian  relaxation  to  solve  a  maximum- 
flow  network-interdiction  problem.  This  problem  arises  when  an  interdictor,  using  limited 
interdiction  resources,  wishes  to  restrict  an  adversary’s  use  of  a  capacitated  network.  In 
particular,  the  adversary,  or  “network  user”  seeks  to  maximize  flow  through  the  network, 
and  the  interdictor  seeks  to  minimize  that  maximum  flow  by  destroying  arcs. 

The  maximum-flow  network- interdiction  problem  (MXFI),  converted  to  an  integer 
program,  is  difficult  to  solve  because  of  an  interdiction-resource  constraint.  Derbes  (1997) 
uses  Lagrangian  relaxation  to  move  this  constraint  into  the  objective  function,  and  is  able 
to  solve  MXFI  efficiently  for  certain  resource  levels  R.  However,  he  solves  for  only  one 
value  of  R  at  a  time.  The  goal  of  this  thesis  is  to  solve  MXFI  efficiently  for  all  reasonable 
values  of  R  in  a  single  procedure. 

The  underlying  reasons  that  motivate  us  to  obtain  solutions  for  all  values  of  R  are: 
First,  R  may  not  be  specified.  For  instance,  at  an  early  stage  of  planning,  we  may  not  know 
how  much  resource  will  be  available  when  the  interdiction  plan  is  to  be  carried  out.  Thus, 
we  would  like  to  provide  the  decision-maker  with  solutions  for  all  potential  values  of  R. 
Secondly,  decision-makers  often  like  to  have  some  input  on  the  amount  of  resource  to  be 
expended.  This  helps  them  review  and  maybe  revise  their  decisions.  Finally,  it  is  not 
much  harder  to  solve  MXFI  for  all  reasonable  values  of  R  than  for  a  single  value. 

Derbes  models  the  interdiction  of  arc  k  as  requiring  a  small,  integer  amount  of 
resource,  rk  >  0.  For  the  sake  of  simplicity,  we  assume  rk  =  1  for  all  interdictable  arcs  k, 

although  arcs  may  also  be  specified  as  non-interdictable.  Derbes  also  considers  a  dynamic 
version  of  MXFI  in  which  flows  require  time  to  move  through  the  network  and  interdicted 
arcs  can  be  repaired  over  time.  The  dynamic  version  of  MXFI  is  beyond  the  scope  of  this 
thesis. 

Bingol  (2001)  investigates  our  version  of  MXFI  and  also  attempts  to  solve  for  all  R 
using  a  Lagrangian-relaxation  procedure.  Essentially,  he  estimates  the  values  of  the 
Lagrangian  functions  for  all  R  simultaneously.  He  assumes  that  these  functions  only  break 
at  points  at  which  the  Lagrangian  multiplier  equals  the  capacity  of  some  arc.  In  carrying 
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out  computations,  he  encounters  some  values  of  R  for  which  his  procedure  cannot  find  an 
exact  solution,  so  he  develops  a  heuristic  to  handle  these  cases  and  detennines  associated 
optimality  gaps.  He  does  not,  however,  optimize  the  Lagrangian-based  lower  bounds. 

In  this  thesis,  we  use  a  procedure  similar  to  Bingol’s  to  find  solutions  for  all  values 
of  R  but  will  optimize  the  Lagrangian  lower  bounds  as  needed.  When  significant 
optimality  gaps  are  identified,  we  apply  a  restricted  branch-and-bound  procedure  to 
improve  the  lower  bounds  and  possibly  identify  better  solutions.  This  procedure  is 
restricted  to  create  enumeration  trees  having  limited  depth. 

A.  BACKGROUND 

There  are  different  types  of  network- interdiction  problems.  For  instance,  Israeli 
(1999)  focuses  on  Maximizing  the  Shortest  Path  (MXSP)  to  slow  a  network  user’s 
movement  between  two  specified  nodes  in  a  network.  The  A-most-vital-arcs  problem  is  a 
special  case  of  MXSP  in  which  the  interdictor  seeks  to  destroy  exactly  k  arcs  to  maximize 
the  shortest  path  length.  Steinrauf  (1991)  designs  a  model  to  isolate  a  targeted  demand 
node  in  a  network.  Morton  (2001)  describes  a  stochastic  network-interdiction  model  to 
help  detect  the  smuggling  of  stolen  nuclear  materials.  However,  in  this  thesis,  we  are  only 
concerned  with  the  detenninistic  maximum-flow  network-interdiction  problem,  MXFI. 

The  scientific  literature  on  network  interdiction  begins  with  Wollmer  (1966)  who 
focuses  on  “Removing  Arcs  from  a  Network,”  where  the  objective  is  to  maximize  the 
reduction  of  the  flow  between  an  origin  and  a  destination  node.  Later,  Wollmer  (1970) 
presents  two  algorithms  for  targeting  strikes  in  an  LOC  (lines-of-communication)  network 
attempting  to  make  the  network  user’s  costs  as  large  as  possible  over  time.  The  LOC 
networks  are  described  by  a  directed  network,  with  arcs  representing  road,  rail,  or 
waterway  segments.  The  costs  in  his  study  include  both  flow  costs  in  dollars,  vehicle 
hours,  manpower  units,  or  any  other  appropriate  measures;  plus  repair  costs  incurred  by 
strikes,  measured  in  similar  units. 

McMasters  and  Mustin  (1970)  devise  an  algorithm  to  detennine  the  optimum 
interdiction  plan  for  minimizing  network  flow  capacity  when  the  minimum  capacity  on  an 
arc  is  positive  and  the  cost  of  interdiction  is  linear  in  the  amount  of  arc  capacity  reduction. 
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Ratiff,  Sicilia  and  Lubore  (1975)  show  how  to  find  a  set  of  arcs,  called  “the  N  most 
vital  links,”  whose  simultaneous  removal  from  a  single-commodity  flow  network  results  in 
the  greatest  decrease  in  the  throughput  capacity  of  the  remaining  system. 

Steinrauf  (1991)  develops  two  mathematical  programs  to  detennine  strategies  to 
interdict  a  network  using  limited  resources.  The  first  model  identifies  a  set  of  arcs  whose 
interdiction  minimizes  the  maximum  flow  through  the  network,  constrained  by  the 
available  resources.  This  is  essentially  MXFI.  The  second  model  identifies  a  set  of  arcs  to 
interdict  that  isolates  a  target  node  and  the  largest  possible  set  of  contiguous  nodes.  The 
latter  model  can  be  used  if  the  exact  location  of  a  target  node  is  uncertain. 

Wood  (1993)  shows  that  the  MXFI  is  NP-complete  even  if  the  interdiction  of  an  arc 
requires  exactly  one  unit  of  resource.  New,  flexible  integer-programming  models  are 
developed  for  the  problem  and  a  number  of  variants  (partial  arc  interdiction,  multiple 
sources  and  sinks,  undirected  networks,  multiple  resources,  multiple  commodities).  Valid 
inequalities  and  a  refonnulation  are  derived  to  tighten  the  LP  relaxations  of  some  of  these 
models. 

Cormican  (1995)  solves  MXFI  using  Bender’s  decomposition,  which  takes 
advantage  of  easy-to-solve  network-flow  subproblems.  Cormican  et  al.  (1998)  solve 
stochastic  versions  of  MXFI  where  arc  capacities  and/or  interdiction  success  may  be 
uncertain.  Stochastic  versions  of  MXFI  are  beyond  the  scope  of  this  thesis,  however. 

Whiteman  (1999)  discusses  techniques  for  planning  the  interdiction  of  a  complex 
infrastructure  network  to  improve  single-strike  effectiveness.  The  motivation  for  his  work 
arises  from  the  huge  cost  of  modern  weapons  and  delivery  platfonns  and  intolerance  for 
casualties,  and  the  consequential  desire  to  minimize  the  number  of  subsequent  attacks.  To 
optimize  target  selection,  he  modifies  the  integer  program  developed  by  Wood  (1993), 
specifying  different  flow  goals. 

Bingol  (2001)  extends  the  earlier  work  by  Derbes  (1997)  on  a  Lagrangian- 
relaxation  technique  to  solve  MXFI  for  all  integer  values  of  total  interdiction  resource,  R, 
in  some  specified  range.  His  basic  procedure  solves  MXFI  exactly  for  most  values  of  R, 
but  “problematic  values”  of  R  arise  with  positive  optimality  gaps. 
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B. 


OUTLINE  OF  THESIS 


This  chapter  has  provided  background  on  network  interdiction  problems  and  has 
surveyed  previous  studies  in  this  field.  The  following  chapter  explains  MXFI  in  detail,  and 
develops  the  Lagrangian-relaxation  technique  for  solving  it,  at  least  approximately.  In 
Chapter  III,  we  first  study  the  algorithms  used  by  Derbes  (1997)  and  Bingol  (2001). 
Subsequently,  we  present  two  new  procedures  to  improve  optimality  gaps  for  problematic 
/?- values  and  sometimes  even  solve  these  problems  exactly.  Chapter  IV  provides  the 
computational  results  for  the  new  procedures.  Finally,  Chapter  V  provides  conclusions  and 
points  out  directions  for  further  research. 
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II.  FUNDAMENTALS  OF  NETWORK  INTERDICTION 


This  chapter  first  defines  necessary  concepts  and  notation  for  our  study.  The 
notation  is  similar  to  that  used  by  Wood  (1993)  and  Bingol  (2001).  The  second  part  of  the 
chapter  starts  with  MXFI,  which  is  a  min-max  model,  and  develops  a  minimizing  integer¬ 
programming  formulation  (MXFI-IP)  and  a  Lagrangian  relaxation  (MXFI-LR)  of  that 
integer  program. 

A.  DEFINITIONS  AND  NOTATION 

G  =  (N,  A)  denotes  a  directed  network  with  node  set  N  and  arc  set  A.  An  arc  k  =  (/, 
j)  starts  from  “tail  node”  i  and  ends  at  “head  node”y.  Each  arc  k  has  a  capacity  Uk  >  0.  The 
network  has  a  source  node  s  and  a  sink  node  t.  For  mathematical  purposes,  we  add  to  A  the 
artificial  arc  a  =  (t,  s)  with  ua=  oo .  The  forward  star  of  node  i,  FS(i),  includes  all  arcs  of 
the  form  (/,  /).  The  reverse  star  of  node  i,  RS(i),  represents  all  arcs  of  the  form  (j,  i ). 

In  this  study,  only  arc  interdiction  is  allowed.  In  the  test  networks,  rk  is  the  amount 
of  interdiction  resource  necessary  to  destroy  arc  k.  For  simplicity,  we  assume  rk  =  1  or  rk 
=  oo;  the  latter  case  simply  means  the  arc  cannot  be  interdicted. 

A  cut  (Ns,  Nt )  partitions  the  node  set  N  into  two  subsets,  Ns  and  Nt  where  5  e  Ns  and 
t  g  Nt.  An  arc  (/,  j)  is  a  forward  arc  of  cut  (Ns,  Nt)  if  i  g  Ns  and  j  g  Nt. 

Ac  represents  the  set  of  forward  arcs  for  cut  (Ns,  Nt).  The  capacity  of  a  cut,  u(Ns, 
Nt),  is  the  sum  of  the  capacities  of  the  forward  arcs  in  that  cut: 

u(Ns,  Nt)=  Yjuk  • 

keAc 

If  u(Ns,  Nt)  is  the  minimum  of  all  cut  capacities,  then  (Ns,Nt)  is  a  minimum  cut. 

We  know  by  the  maximum-flow  minimum-cut  theorem  that  the  capacity  of  the  minimum 
cut  equals  the  maximum  s-t  flow  in  that  network;  the  maximum  flow  model  will  be 
described  shortly. 
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B. 


NETWORK  INTERDICTION  MODEL 


Here,  we  define  the  basic  MXFI,  which  is  then  converted  into  MXFI-IP. 
Lagrangian  relaxation  for  MXFI-IP  is  presented  subsequently. 

1.  Maximum  Flow  Network  Interdiction  Problem  (MXFI) 

In  MXFI,  an  interdictor,  using  limited  interdiction  resources,  destroys  arcs  of  the 
enemy’s  network  G=  (N,  A)  and  makes  them  useless.  The  objective  is  to  minimize  the 
maximum  amount  of  flow  that  the  enemy,  i.e.,  the  network  user,  can  transport  through  the 
network.  We  define  MXFI  as  a  function  of  total  interdiction  resource  R : 


MXFI  (R): 

Indices: 

i,jeN 
s  e  N 
t  e  N 

k  =  ( i,j )  e  A 
a  =  (t,  s )  eA 

Data: 

Uk 

rk 

R 

Decision  Variables: 

yk 

xk 


nodes  of  network  G={N,A ) 
source  node 
sink  node 
network  arcs 

artificial  return  arc  included  in  A 

nominal,  uninterdicted  capacity  of  arc  k 

amount  of  resource  necessary  to  interdict  arc  k  (ra  =  cc) 

total  amount  of  available  resource 

amount  of  flow  on  arc  k  (network  user) 

1  if  interdiction  of  arc  k  is  attempted  (interdictor) 

0  otherwise 
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Formulation: 


z(R )  =  min  max  ya 

xeX  y 


Zw  "  Z-!h-=0 

V/  e  N 

(1) 

keFS(i)  keRS(i) 

\/k  e  A 

(2) 

where  X  =  i  x  e  {0,  l}^ 

Yjrkxk  ^  R>xa  =  °|- 

(3) 

l 

keA  J 

In  this  model,  constraints  (2)  force  capacities  down  to  zero  for  all  interdicted  arcs. 
For  fixed  x,  the  inner  maximizing  problem  is  a  standard  maximum-flow  model,  which 
maximizes  the  flow  from  a  source  node  to  a  sink  node  subject  to  flow-balance  constraints 
(1)  and  arc  capacities  (2).  MXFI(7?)  can  also  be  reformulated  as  MXFIR(R)  (Cormican  et 
al.  1998): 


MXFIR(R): 


z(R)  =  min  max  y  -^xkyk 

xeX  y 


keA 


s-t.  -  2>a-=0 

keFS(i)  keRS(i) 


Vi  eN 


0  ^  yk  ^  uk 


\/k  e  A 


where  X  =  j  x  e  {0,  l}^  ^ rkxk  <R,xa=  0 


keA 


2.  An  Integer  Program  for  MXFI  (MXFI-IP) 

Taking  the  dual  of  the  inner  maximization  of  MXFIR  gives  us  a  linear  minimization 
model  (which  is  derived  by  Wood  1993  directly  from  MXFI): 
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MXFI-IPCR) 
Indices:  As  in  MXFI 


Data:  As  in  MXFI 
Decision  Variables: 

a,-  1  if  /  e  Nt  and  0  if  i  e  Ns,  for  a  model-identified  cut  ( NsJVt ) 

jik  if  A:  is  a  forward  arc  of  cut  (Ns,Nt)  and  not  interdicted,  0  otherwise 

%k  1  if  interdiction  of  arc  k  is  interdicted,  0  otherwise  (interdictor) 

Formulation: 

z(R)  =  min  £  ukpk 
“•p-* 


s.t. 

a,  -  a}  +xk+pk  >0 

\fk  =  (ij)  e  A  -{a} 

(4) 

a,  ^ 

for  a  =  ( t ,  s) 

(5) 

Ysrkx*-R 

(6) 

keA 


a.i  e  {0,1}  V/  e  N 
xk,  fik  e  {0,1}  \/k  e  A 
as  =  0,  at  =  1 
Pa  =  %a  =  0 

3.  Lagrangian  Relaxation  for  MXFI  (MXFI-LR) 

Lagrangian  relaxation  is  a  general  solution  strategy  for  solving  mathematical 
programs  that  pennits  us  to  decompose  problems  to  exploit  their  structure.  This  solution 
approach  leads  to  bounds  on  the  optimal  objective  function  value  and,  frequently,  to  good, 
though  not  necessarily  optimal  solutions  (Ahuja,  Magnanti  and  Orlin  1993,  p.  601). 
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MXFI-IP  can  be  difficult  to  solve  because  of  constraint  (6),  so  we  transfer  this 
constraint  into  the  objection  function  following  Derbes  (1997).  After  this  modification,  the 
new  model,  MXFI-LR,  can  be  easily  solved,  although  it  may  not  always  yield  an  optimal 
solution.  The  model  is: 


MXFI-LR  (A,  R): 


z(A,R )  =  min 

a,p,x 


YjlhPk  +  M£rkxk-R) 

keA  keA 


s.t.  a i  -  a  j  +  xk  +  Pk  >  0 


\/k  =  (ij)  e  A  -  {a} 


a,  -  as  +  xa  +  j3a  >1  for  a  =  (t,  s ) 
a,-  e  {0,1}  V/  e  N 
Xk,  pk  e  {0,1}  \/k  e  A 
as  =  0,  a,  =  1 
Pa  =  %a  =  0 


(V) 


The  term  A,(^rk x k  -  R)  represents  a  penalty  for  exceeding  or  a  reward  for  not  consuming 

ksA 

the  resource  budget  R.  The  parameter  A  can  be  interpreted  as  the  dual  cost  of  interdicting 
an  arc. 

The  linear-programming  relaxation  of  MXFIR  has  integer  extreme  points,  so  we 
can  take  the  dual  of  that  relaxation  to  obtain  MXFI-LRD  after  some  simplifications 
(Derbes  1997): 


MXFI-LRD  (A  ,R): 

z(A,  R)  =  max  ya  -  AR 

y 

s-t.  'Ey*  -  Z>k=0  V/eA 

keFS(i)  k<=RS(i) 
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0  <  yk  <  min  {uk ,  Ark  }  \/k  e  A 


Of  course,  z(A,R )  is  only  a  lower  bound  on  z(R),  but  we  can  use  this  lower  bound, 
and  interdiction  plans  obtained  during  its  solution  to  solve  MXFI(f?),  at  least 
approximately.  Procedures  to  do  this  are  presented  in  the  next  chapter. 
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III.  NETWORK  INTERDICTION  BY  L AGRAN GIAN  RELAXATION 

AND  BRANCH-AND-BOUND 


This  chapter  explains  our  algorithm  for  solving  MXFI.  The  algorithm  uses  MXFI- 
LRD  as  Bingol  (2001)  does,  but  adds  two  new  procedures  to  decrease  optimality  gaps.  The 
first  procedure  finds  optimal  Lagrangian  multipliers,  which  Bingol’s  approach  does  not  do. 
The  second  procedure  implements  a  restricted  version  of  branch-and-bound.  “LRD”  and 
“MXF1-LRD”  are  used  interchangeably  below. 

A.  PRELIMINARIES 

Derbes  (1997)  shows  that  Lagrangian  relaxation  can  fail  to  solve  MXFI(R)  when 
capacities  are  not  unique,  but  may  solve  the  problem  correctly  if  the  capacities  are 
perturbed  slightly  so  as  to  be  unique.  We  use  the  procedure  suggested  by  Bingol  (2001)  to 
guarantee  unique  capacities: 

Procedure  Perturb(u): 

{ 

Order  uk  so  that  ux  <  u2<  —  <  u\a\  ; 

For  ( k  and  m  such  that  Uk-\<  Uk  Uk+ 1  Uk+m<  Uk+m+l)i 

Uk<—  104  Uk, 

uk+l  <—  10  Uk  +  1 


uk+m<~  uk  +  m  -  1; 

} 

Return  u; 


So,  suppose  that  u  =  (1,  1,  3,  3,  3).  After  running  this  procedure,  the  network  has 
unique  arc  capacities,  u=  (10000,  10001,  30000,  30001,  30002).  These  small  perturbations 
are  unlikely  to  change  the  optimal  solutions  to  MXFI.  (A  scaling  factor  larger  than  104 
would  have  to  be  used  if  the  network  were  very  large  and  there  were  many  arcs  with 
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identical  original  capacities.)  The  use  of  this  procedure  is  assumed  and  not  shown 
explicitly  below. 

B.  THE  LAGRANGIAN-RELAXATION  MODEL 

The  function  z(A,R),  from  MXFI-LRD(T,R) ,  is  a  piecewise-linear  concave 
function  in  A  for  fixed  R.  Since  LRD(T,R)  relaxes  MXFI-IP(R),  z(A,R)<z(R)  for  all 
A ,  where  z(R)  is  the  optimal  solution  for  MXFI-IP(R).  Naturally,  we  wish  to  maximize 
z(A,R)  overT . 

Given  a  fixed  A ,  we  can  solve,  or  attempt  to  solve,  MXFI(R)  and  obtain  the 
interdiction  set  Af(A)  ,  for  an  unspecified  R,  through  LRD  (A,  R)  as  follows  (Bingol  2001): 

1 .  Modify  G  =  ( N ,  A)  such  that  arc  capacities  become  u[  =  min)^ ,  Ark }  \/k  e  A, 
and  use  a  max-flow  algorithm  to  find  a  minimum-capacity  cut  Ac , 

2.  Define  the  interdiction  set  AZ(A )  =  {k  e  Ac  \  uk  >  Ark  }. 

The  solution  found  through  LRD  (A,  R)  identifies  a  minimum  cut  ( Ns ,  Nt )  that 
corresponds  to  a  feasible  or  infeasible  solution  x  to  MXFI(R),  where  xk  =1  \/k  e  Af(A) , 
and  xk  =  0  otherwise.  If  x  satisifies  the  budget  constraint  ^  rkxk  <  R ,  then  it  is  a  feasible 

keA 

solution.  By  solving  LRD(/t,f?)  for  several  values  of  A,  we  may  be  able  to  find  x, 
with  '^>\xk  =  R  .  (Intuitively,  this  seems  more  likely  to  occur  at  A  =  uk  which  is  why 

k&A 

Bingol  only  evaluates  such  values  of  A.)  In  this  case,  MXFIff?)  is  solved  exactly.  If  not, 
we  will  establish  a  non-zero  optimality  gap. 

To  determine  optimality  gaps,  we  need  upper  and  lower  bounds  on  z(R) .  For  a 
feasible  interdiction  set  Af(A)  ,  the  corresponding  maximum  flow  is  an  upper  bound  (UB); 
z(A,  R)  from  LRD  is  a  lower  bound  (LB)  for  z(R)  for  any  A;  the  best  LB  from  the 
Lagrangian  relaxation  is  found  by  solving  max  z(A,R) . 

A 

Derbes  (1997)  uses  a  binary  search  on  the  interval  of  uncertainty  for  A  to  maximize 
z(A,R ) .  His  algorithm,  when  specialized  for  our  version  of  MXFI(R),  is: 
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Algorithm  D(i?): 


Input: 


G=  (N,  A)  directed  network  with  source  seN  and  sink  teN;  arcs  indexed  by  k, 
uk  unique,  scaled-up  integer  arc  capacities  uk  >  0 ,  Vk  e  A, 


n 

R 

Output:  LB 


rk=  1  for  all  interdictable  arcs  k  e  A;  rk  =  oo  otherwise, 

total  interdiction  resource  available. 

lower  bound  for  the  network  interdiction  problem  given  R ; 


{ 


} 


Initialize:  LB  < - oo  :  A  <—  u  +  1:  A  <—u  —1 :  A  <— A  : 

while  ( Amax  -  Amin  >  1  )  { 

u'k  <—  min \uk,Ark}  for  all  k  e  A ;  /*  adjust  arc  capacities  */ 

/*  maxFlow(  )  below  returns  a  min  cut  Ac  and  max  flow  value  /  */ 
( Ac,f )  <—  maxFlow  (G,  s,  t ,  u' ); 

{k  e  ylc  I  «[=  /l  }; 

LB  <—  max{LB,f  -  AR) ; 

if  (|  Aj  |  =  f?)  {  break;  }  /*  from  while  loop  */ 

if  (I4I<*) 
if  (1 4 1 > R) 


^  +^max)/2  i 

} 

print  (“  LB  for  f?  =”,  /?,  “is”,  Z5); 


In  Algorithm  D(f?),  the  interval  of  uncertainty  initially  has  left  endpoint 
4m  =  - 1  and  right  endpoint  Anmx  =  wmax  + 1 .  For  A  =  (Aimx  +AmJ/2,  a  max-flow 

problem  is  solved  in  G  =  (N,A )  with  u'k  =  max {uk,Ark}  and  A,  is  determined.  The 
endpoints  of  the  interval  are  adjusted  by  detennining  the  sign  of  the  slope  of  z(A,R)  which 
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is  sgn (^rkxk  -  R)  :  If  |  A{  \  <  R  ,  2  becomes  the  right  endpoint,  2max ,  of  the  interval.  If 

k 

\Aj\>  R,  then  X  becomes  the  left  endpoint,  Xmin ,  of  the  interval.  If  \Ak\  =  R ,  then 
MXFI(f?)  has  been  solved  optimally.  If  an  optimal  solution  to  MXFI(f?)  is  not  found,  this 
process  goes  on  until  X  is  found  such  that  2  .  <  X  <  2m„  and  2v  -  X  .  <  1 .  This  X  is 

close  enough  to  X* ,  because  the  arc  capacities  have  been  scaled  up  by  a  large 
multiplicative  factor.  (See  Section  A  of  this  chapter.) 

With  this  procedure  we  can  obtain  a  solution  for  only  a  single  value  of  R.  But,  we 
are  interested  in  solving  MXFI(R)  for  all  “reasonable”  values  of  R,  which  means 
f?e[0,|  Amjn  |],  where  Amin  is  a  minimum-cardinality  cut  consisting  of  only  interdictable 

arcs.  (We  assume  such  a  cut  exists.)  In  other  words,  \Amm\  is  the  minimum  resource  level 
that  guarantees  a  feasible  interdiction  plan  can  force  the  network  flow  to  zero. 

Perfonning  a  binary  search  for  each  R  e  [0,  |  Amin  |]  seems  inefficient.  This  is 
especially  true  because  z(X,R)  and  z(2,f?')  are  closely  related: 

z(X,R)  =  z(X,R')  +  X(R'-R). 

Furthermore,  z(X,R)  =  z(X,0)-XR  where  z(2,0)  is  simply  the  maximum  flow  f(X)  in 
G  =  ( N,A )  with  arc  capacities  u[  =  max {uk,Xrk}  \/k  e  A. 

Bingol  (2001)  attempts  to  solve  MXFI(f?)  for  all  (reasonable)  values  of  R  by 
approximately  solving  LRD  (X,  R) .  He  considers  X  in  the  range  [0,wmax  ]  and  assumes  the 

slope  of  z(X,R)  =  / (2)  -  XR  can  change  only  at  X  =  uk .  That  is,  he  assumes  that  Af(X ) 
changes  only  in  going  from  X  =  uk-£  to  2  =  uk  +  s  for  certain  uk  and  some  small  s  >  0  . 

With  this  assumption,  he  can  solve  many  instances  of  MXFI(f?)  with  acceptable  accuracy 
by  only  evaluating  z(X,R)  at  such  points.  Bingol’s  algorithm  is  given  next.  The  routine 
ProblematicRn()  initially  represents  Bingol’s  procedure  ProblematicRl()  for  handling 
problematic  R-values.  That  will  be  replaced  subsequently  with  our  routine, 
ProblematicR2(),  to  improve  results. 
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Algorithm  A:  Attempts  to  solve  MXFI(R)  for  each  value  of R  e  [0,  Amin  |]  by  computing 
z(A,R )  for  A  =  uvu2,...,ulAl  and  determining  a  feasible  interdiction  set  A\R) 


Input:  G=  (N,  A)  directed  network  with  source  ssN  and  sink  teN;  arcs  indexed  by  k 


n 

Output:  x(R) 

A}*> 

z(R) 

UB(R) 
LB(R ) 


unique,  scaled  up  integer  arc  capacities  uk  >  0  ,  Vk  e  A, 
r/c=  1  for  all  interdictable  arcs  k  e  A;  rk  =  go  otherwise, 
optimal  or  suboptimal  solution  for  MXFI(A), 
interdiction  set  V  R  e  [0,  |Amin|]  that  corresponds  to  x(R) , 

for  non-problematic  R:  optimal  objective  function  value  to  MXFI(A), 
for  problematic  R:  feasible  objective  value, 
for  problematic  R:  lower  bound  on  optimal  objective  value, 
absGap(R )  for  problematic  R:  absolute  optimality  gap  UB(R)  -  LB(R), 
relGap(R)  for  problematic  R:  relative  optimality  gap  \00xabsGap(R)/LB(R). 


/*  define  and  order  the  candidate  set  of  A  values  */ 

Reorder  arcs  so  that  u,  <u2<  ■■■  <  u,A, ; 

Initialize:  R  <—  |A|+1;  R'  <—  |A|+1;  /  <— 0;  /' <— 0;  A'  <—  0;  /  <—  1; 

while  ( /  <  \A  \  and  R  -*■  0)  { 

A  ^ —  Uj  \ 

u[  <—  min{ uk,Ark)  for  all  k  e  A ;  /*  adjust  arc  capacities  */ 

( Ac,f )  <—  maxFlow  (G,  s,  t,  u' ); 

Aj <—  {k  s  Ac  |  u'k=  A  }; 

xk  <—  1  Vk  e  Aji  xk  < —  0  Vk  e  A\  Ar; 

R  <r-  |  A/\; 

z  <—  f-  AR ; 

/*  ProblematicRl()  is  called  by  Bingol;  we  will  call  ProblematicR2()  */ 
if (R'  - R>  1  ){  ProblematicRn(R',R,/',/,/l',2,ylj);  } 

if  (R  *  R')  { 
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print  (“  The  optimal  solution  to  MXFI(i?)  for  R  =  ”,  R,  “is  ”); 
print  (x,A;,  z)  ; 

} 

/*  In  the  following  case,  we  know  that  MXFI(f?-l)  will  be  solved  by 
evaluating  z(A  +  s,R- 1)  for  a  small  s  >  0  */ 
if  (  uk,  =  A  for  some  k '  e  Aj){ 

Ai<—  Aj-  {A:'}; 

z  <—  z  +  uk,  ; 

xkf  <—  0; 

print  (“  The  optimal  solution  to  MXFI(f?)  for  R  =  ”,  R,  “is  ”); 
print  (x ,  A /,  z)  ; 

} 

if  (R'*R)  { 

R'<-R;  A',  <-  A,  z'^z; 

} 

/  <-  /  + 1 ; 

} 

} 


Algorithm  A  orders  arcs  by  increasing  capacity  ux  <  u2  <  ■  ■  ■  <  u  A  and  starts  with 
A=uv  Because  all  arc  capacities  are  A,  AI  =  Ac,  which  corresponds  to  a  post¬ 
interdiction  flow  of  zero.  This  effectively  defines  the  solution  x(r;  )  to  MXFI(rj)  ,  where  rx 
is  the  number  of  arcs  in  a  minimum-cardinality  cut  consisting  of  only  interdictable  arcs.  If 
uk,  =  A  for  some  k'  e  A/,  we  also  find  a  solution  to  M X F I (  / ,  -1).  If  not,  the  next  iteration 

sets  A  to  u2  and  attempts  to  find  a  solution  x(r2)  to  M  X  F I  (>2 )  such  that  r2  <  rx .  If 
r2  =  rx  - 1 ,  the  algorithm  continues  without  need  of  ProblematicRn(),  and  repeats  for 
increasing  values  of  A  =  uk . 
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When  we  compute  R  =  |  At  \ ,  we  know  that  we  solve  max  z(X,\Af  |)  =  z(X*r,R)  , 

X 

i.e.,  we  maximize  the  Lagrangian  function  z(X,R)  for  R  =  \Ar\.  Furthermore,  since 
^rkxk  =R,  z(X*r,R)  =  z(R) ;  therefore  we  have  solved  MXFI(R).  However,  difficulties 

ksA 

arise  when  the  algorithm  finds  z(X*r,R)  =  z(R)  and  z(X*r,,R')  =  z(R')  for  R' >  R  + 1,  but 
not  for  R  +  l,R  +  2,...,R'  We  call  these  “missing”  or  “problematic  i?-values.”  For 
problematic  R- values,  Bingol  (2001)  derives  a  heuristic,  ProblematicRl()  to  find  a  feasible 
solution  and  lower  bound  LB  on  z{X*r,R) .  The  feasible  solution  gives  an  upper  bound  UB 
for  z(R) ,  so  he  can  compute  an  absolute,  or  relative  optimality  gap  which  are  defined  as 
absGap=UB~LB  and  relGap=  100  x  absGap/LB,  respectively. 


Procedure  ProblematicRl  (i?',R,/',/,2',2,,^) 

Input:  R  current  R-value  correctly  solved, 

R'  previous  R-value  correctly  solved  (R'  >  R  +  1), 

X  value  of  parameter  that  maximized  z{X,R)  to  solve  MXFI(R), 

X'  value  of  parameter  that  maximized  z(X,R')  to  solve  MXFI(i?'), 

/  maximum  flow^^-), 

/'  maximum  flow  / ( X' )  , 

A'j  interdiction  set  for  MXFI(  R' ). 


Output:  \(R) 

4S) 

UB(R ) 
LB(R) 
absGap{R) 
relGap(R) 


feasible  (possibly  optimal)  solution  for  MXFI(R), 
interdiction  set  that  corresponds  to  x(i?) , 
feasible  objective  value  for  x(R) , 
lower  bound  on  optimal  objective  value  z(R), 
absolute  optimality  gap  UB(R )  -  LB(R )  for  x(R) , 
relative  optimality  gap  100 xabsGap(R)/LB(R)  )  for  x(R) . 


Initialize:  UB  <—  /'  -  XR' ;  Ak  <—  A' ; 
for  (  r  =  R’~  1  down  to  R  +  1  )  { 
kmm  <-  argmin  uk ; 

keA} 

LB  <— max{/-  X r,f-  X'r  }; 

UB  <—  UB  +  uk  ; 

^min 

absGap  <—  UB  -  LB; 
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relGap  <—  100 xabsGap/LB', 

A  A  ~  ^min  ’ 

xk  <—  1  \/ke  Aj\  xk  <—  0  \/k  e  A\Af; 

print  (“  For  problematic  R  =  ”,  r,  “the  approximate  soln.  to  MXFI(f?)  is:  ”); 
print  (x ,  Aj  ,UB ,  LB,  absGap,  relGap  ); 

} 

} 


Algorithm  A  requires  the  solution  of  a  sequence  of  maximum- flow  problems.  We 
use  the  shortest  augmenting  path  max-flow  procedure,  maxFlow(),  (Edmonds  and  Karp 
1972),  which  finds  the  maximum  flow/  and  a  minimum-capacity  cut  Ac-  Since  the  flow 
on  each  arc  is  non-decreasing,  this  procedure  uses  the  previous  flow  level  from  the  last 
iteration,  instead  of  starting  from  scratch.  This  is  important  for  fast  run-times. 

C.  FINDING  AN  OPTIMAL  A  FOR  PROBLEMATIC  R 

The  procedure  ProblematicRl()  above  does  not  find  the  best  LB,  i.e.,  it  does  not 
maximize  z(A,R)  for  problematic  /Avalucs  and  therefore  does  not  establish  the  best 
achievable  optimality  gap.  We  can  attempt  to  find  the  best  lower  bound  by  exploring  the 
region  by  defining  a  finer  grid  for  A  .  But  this  is  computationally  expensive  and 

requires  solving  many  max-flow  problems.  At  this  point,  Proposition  1  expedites  the 
process. 

Proposition  1  (Derbes  1997):  The  Lagrangian  function  z(A,R)  has  a  slope  of 
|  Aj(A)  |  -R  ,  where  |  Af(A)  |  is  the  number  of  arcs  interdicted  using  A . 
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a  9 


e  10 


Figure  1.  NETSample.  Letters  and  numbers  represent  arc  labels  and  capacities,  respectively. 
Unlabeled  arcs  cannot  be  interdicted.  All  labeled  arcs  have  rk  =  1 . 


Figure  1  shows  a  sample  network,  NETSample,  which  we  will  use  to  explain 
Proposition  1.  We  will  compute  z(2,  l)and  look  at  the  function  in  detail: 


Figure  2.  Lagrangian  function  for  z(A,  1)  for  NETSample.  The  first  number  in  parentheses  is  the 
value  of  parameter  A  and  the  second  number  represents  the  value  of  the  function. 

If  we  compute  the  slope  mBc  of  line  BC  in  Figure  2,  we  see  that  mBC  =  \  Aj  (2)  |  -R 
=  3-1=2.  Likewise,  mCD  =  \  Aj(A)  |  -R  =  1-1  =  0 ,  which  also  gives  us  the  optimal 
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solution  for  R  =  \ .  (When  the  slope  is  zero,  \AI(A)\  =  R  =  \,  and  we  obtain  an  exact 
solution  for  MXFI(l).)  These  results  are  in  accordance  with  Proposition  1. 

For  NETSample,  if  we  run  Algorithm  A,  we  run  into  a  missing  value  of  R.  We  find 
a  solution  for  R  =  3  at  ^  =  3  and  R  =  1  at  \  =  7 ,  but  not  for  R  =  2  .  The  solutions  for 

R  =  3  and  R  =  1  are  A(f3)  =  {a,  b,  c}  and  A31'  =  {e}.  Bingol’s  heuristic  finds  the  UB  and 
LB  for  R  =  2  ,  with  the  following  steps: 

i.  42)  =  A^-iknun)  =  {a,b,c}-{c}  =  {a,b}, 

ii.  UB  =  z(  3)  +  umin  =  1  +  7  =  8, 

iii.  LB  =  max {/(/Ij ) - 2/lj , /(/l3 ) - 2/l3 }  =  max  {19-14,  10-6  }  =  5, 

iv.  absGap  =  8-5  =  3, 

v.  relGap  =100  x  3/5  =  60%. 

The  relative  gap,  relGap,  is  quite  high.  To  reduce  this  optimality  gap,  we  exploit 
Proposition  1  and  trace  the  function  z(A,R)  more  closely. 
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Figure  3.  Lagrangian  Function  of  R  =  2  for  NETSample.  The  solid  line  denotes  that  the  value  of  the 
function  is  computed  for  those  A  values.  The  dotted  line  shows  how  procedure 
solveExactQ  computes  the  optimal  A  . 


ProblematicRl()  finds  LB  =  5  for  R  =  2,  which  is  the  value  of  the  function  at  point 
E.  But  with  Proposition  1,  we  can  look  for  /C  over  the  region  [3,  7]  to  improve  LB  and 
possibly  solve  MXFI.  If  Af(AD)  =  2  (it  is  not  in  this  case),  then  MXFI(2)  is  solved 
exactly.  Otherwise,  we  obtain  an  improved  LB.  To  accomplish  this,  we  do  the  following: 

i.  /*  Find  the  slopes  of  associated  lines  */ 

mCD  =  |  Aj(3)  |  —R  =  3-2  =  1,  mDE  =  \  Af(7)  |  - R  =  1-2  =  -  1, 

ii.  /*  A*  is  the  optimal  parameter  value  for  R  =2  */ 

z(  A*2,2)  =4  +  (1)(  A*2-  3),  z  ( A*, 2)  =  5  +  ( -1 )( X-  7), 

iii.  /*  Compute  A*2  from  the  equations  above  */ 

A*  =5.5  and  z  (A*, 2)  =  6.5, 

iv.  /*  Check  for  optimality  */ 

Aj (5.5)  ^  2  ,  /*  MXFI(2)  is  not  solved  exactly  */ 
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v.  /*  Compute  the  optimality  gap  with  the  new  LB  */ 
absGap=  UB-z(A^,2)  =  8-6.5  =  1.5,  and 
relGap  =  100 xabsGap / z(/ t,*,2)  =  100  x  1.5  /  6.5  =  23.1%. 


The  relative  gap  has  decreased  from  60%  to  23.1%.  The  best  Lagrangian  multiplier 
is  A*2  =5.5.  Because  Algorithm  A  computes  z{A,R)  at  A  =  uk\/k  e  A,  it  would  never 

discover  this  LB.  This  example  has  only  one  missing  R-value.  However,  there  can  be 
contiguous  problematic  R-\ alues.  The  following  procedure  will  also  handle  these  cases. 


Procedure  ProblematicR2( R\ R , /', /,  A', A,  A)): 


Input:  R 

R' 

A 

A' 

f 

f 


current  R-value  correctly  solved  , 
previous  R-value  correctly  solved  (R'  >  R  +  1), 
value  of  parameter  at  which  R  is  solved, 
value  of  parameter  at  which  R'  is  solved, 
maximum  flow  /(  A  ), 
maximum  flow/(  A' ). 


Output:  \(R) 

4R) 


optimal  or  suboptimal  solution  for  MXFI(f?), 
interdiction  set  V  R  e  [0,  |Amin|]  that  corresponds  to  x(R) , 


z(R)  for  non-problematic  R:  optimal  objective  function  value  to  MXFI(R), 

UB(R)  for  problematic  R:  feasible  objective  value, 

LB(R)  for  problematic  R:  lower  bound  on  optimal  objective  value, 
absGap(R )  for  problematic  R:  absolute  optimality  gap  UB(R )  -  LB(R), 
relGap(R)  for  problematic  R:  relative  optimality  gap  100 xabsGap(R)/LB(R). 

For  (r  =R'  -  1  to  R+l)  { 

point  <—  A; prevPoint  <—  A' ;  Rprev  <- R'  ;  Rcur  R' ; 


if(r  =  R  +  l){Rci„.^R;} 
while  {point  L-  prevPoint)  { 
prevPoint  <—  point ; 

point  <—  (/  -  /'  -  AR  +  A’R')I{R'  -R);  /*  see  the  explanation  below  */ 
uk  <—  min{MA. ,  {point)  rk  }  for  all  k  e  A; 
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{Ac,  flow)  <—  maxFlow(G,  s,  t,  u' ); 

A  i  ^ { k  (£  A  q  a  *  =  2  }  \ 

xk  <—  1  \/k  e  Aj ;  xk  <—  0  \/k  e  A\Ak; 

/*  If  the  exact  solution  is  found,  then  it  prints  the  results  */ 
if  (|4|  =  r)  { 

print  (“  The  optimal  solution  to  MXFI(R)  for  R  =  ”,  R ,  “is  ”); 
print  (x,A/,z)  ; 

prevPoint  <—  point ;  /*get  out  of  while  loop  */ 

4»r  r  ’ 

^  (  R  prev  ~  Rcur  =  1  )  {  R  prev  5  } 

} 

if  ( |  Aj  |  >  r  )  {  2'  <—  point ;  R'  <—  |  4  I  5  /'  <—  flow,  A\  <—  4  i  } 
if  ( |  4  I  <  r  )  {  ^  <—  point ;  f?  <—  1 4  U  /  fl°w'f 

} 

4  <—  point ;  /*  The  best  Lagrangian  multiplier  is  found  */ 

if(^-4,,>D{ 

/*  With  revised  values  of  2  and  /l' ,  we  can  reuse  ProblematicRl()  to 
find  a  feasible  solution  and  optimality  gap  for  missing  values  of  R  */ 
ProblematicRl  ( Rprev ,  Rcur ,  /',  /,  2',  2, 4V"  ’ ); 

/*  Uncomment  the  next  statement  to  invoke  Branch()  .  Notice  that  UB(R) 
comes  from  ProblematicRl ()  */ 

/*  Branch(  ,  Rcur,  A*+l  -  s ,  A^ ,  UB(R));  */ 

D  > _  n 

prev  vcur  5 

} 

} 

} 

In  our  heuristic  above,  we  derive  the  formula  point  =  (f  -  f  -  AR  + A' R')  l(R'  -  R) 
in  order  to  find  the  intersection  with  the  following  steps: 
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1.  /*  Find  z(point,r )  using  the  previous  solution’s  parameters  ( R',A',f '). 
Notice  that  r  is  the  missing  value  and  slope  is  R'  -r  */ 

z(point,  r)  =  (/'  -  A'r)  +  ( R '  -  r)(point  -  A') , 

2.  /*  Find  z(point,r )  using  the  current  solution’s  parameters  ( R,A,f).  This 
time,  slope  is  r  -  7?  */ 

z{point,  r)  =  (/  -  Tr)  +  (r  -  R){A  -  point ) , 

3.  /*  Using  these  two  equations,  solve  for  “ point' ’  */ 
point  =  {f-f-AR  +  A'R')/(R'  -R). 

When  Algorithm  A  runs  into  one  or  more  problematic  R -values,  ProblematicR2()  is 
called  to  find  the  best  Lagrangian  multiplier  for  each  missing  value  of  R.  This  heuristic 
seeks  T*  over  an  interval  of  uncertainty  by  performing  a  binary  search  based  on  the  sign  of 

slope  ^rkxk  -  R  ■  First,  it  finds  the  intersection,  point,  using  the  fonnula  derived  earlier. 
If  Aj  {point)  =  r  ,  then  MXFI(r)  is  solved  exactly.  Otherwise,  it  updates  the  endpoints  as  in 
Derbes’s  procedure  (1997).  It  keeps  updating  until  it  converges  to  A*R  and  possibly  finds 
an  exact  solution.  If  there  are  multiple  values  of  R  e  [R  - 1  ,RCW.  + 1] ,  ProblematicRl()  is 

invoked  to  improve  the  results  with  the  best  Lagrangian  multipliers  for  each  of  these 
problematic  R -values.  For  further  improvements,  we  can  a  invoke  Branch()  which  applies 
a  restricted  branch-and-bound  procedure. 

If  there  is  only  one  missing  /?- value,  one  of  two  cases  occurs,  and  it  only  takes  one 
evaluation  of  z(A,R)  to  figure  this  out: 

1 .  If  there  is  only  a  single  breakpoint,  then  we  must  compute  it  exactly  when 
we  compute  the  intersection  of  the  two  lines  extending  upwards  from  the  coordinates 
{Ar,,z{Ar.,R'))  and  ( Ar,z{Ar,R ))  in  {A  ,  z)-space.  Let  that  intersection  point  be  (A"  ,z"). 

The  single  breakpoint  is  verified  by  computing  z(A",R )  and  finding  z"  =  z(A",R )  . 

2.  If  there  are  two  breakpoints  in  the  interval  {AR.,AR),  we  compute  (A"  ,z"), 
evaluate  z{A",R )  and  find  that  z{A",R )  <  z" .  Consequently,  the  slope  of  the  z{A,R )  curve 
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must  lie  strictly  between  one  and  negative  one.  Since  the  slope  can  only  change  by  integer 
amounts,  we  have  found  a  portion  of  the  curve  where  the  slope  is  zero  and  will, 
consequently,  have  solved  MXFI(R)  optimally. 

If  there  are  sequential  problematic  R -values,  there  are  three  possibilities: 

1.  There  may  not  be  an  exact  solution  for  any  of  the  problematic  //-values.  In 
this  case,  the  procedure  finds  A*  and  then  ProblematicRl()  is  called  with  updated 
parameters  to  find  the  best  LB.  In  fact,  A  and  A'  computed  earlier  by  Algorithm  A  are  at 
some  uk .  However,  when  computed  by  ProblematicR2(),  either  A  or  A'  becomes  the  best 
Lagrangian  multiplier,  which  helps  improve  LB  for  these  problematic  //-values.  To  invoke 
Branch(),  we  use  the  following  idea:  When  there  are  more  than  one  problematic  R -values  at 
one  breakpoint,  the  best  Lagrangian  multiplier  is  the  same  for  all  these  values.  In  other 
words,  for  R  e  [R' -  1,R  + 1] ,  A*+1  =  A*R+2  =■■■  =  A*  , .  Therefore,  we  use  only  one  of  them 

with  a  slight  perturbation,  A*+1  -  s ,  to  invoke  Branch(),  where  |  Af  (A*+l  -£•)!=  R  =  R' . 

2.  The  interval  of  uncertainty  might  yield  an  exact  solution  for  some  of  the 
problematic  //-values.  If  it  does,  then  the  procedure  prints  out  the  results  of  those  values, 
and  R  and  Rcur  establish  new  endpoints  for  the  problematic  //-values.  If  R  -  Rcur  >  1, 

ProblematicRl()  and/or  Branch()  can  be  invoked  to  find  a  feasible  solution  and  optimality 
gap  for  problematic  R  e  [Rcur  + 1,  R  - 1]  as  explained  in  the  previous  paragraph.  Since 

Ai/R,mv>  and/or  A\Rcur)  might  change  for  the  problematic  R  e  [Rrur  + 1,  Rprev  - 1] ,  not  only  LB 
but  also  UB  are  computed  by  ProblematicRl(). 

3.  The  procedure  might  solve  MXFI(R)  exactly  for  all  missing  //-values  in  the 
sequence.  In  our  tests  in  Chapter  IV,  we  observe  this  in  two  instances  out  of  10  test 
problems. 

D.  RESTRICTED  BRANCH-AND-BOUND 

Using  procedure  ProblematicR2()  to  solve  MXFI,  we  may  be  able  to  reduce  the 
optimality  gap  by  a  substantial  amount  or  solve  MXFI  exactly,  for  certain  values  of  R. 
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However,  there  might  be  some  problematic  /Lvalues  for  which  the  optimality  gap  is  still 
large.  For  these  R -values,  we  can  resort  to  branch-and-bound. 

The  usual  idea  of  the  branch-and-bound  in  the  setting  of  integer  programming  is  the 
following:  Assume  that  we  have  a  feasible,  binary  integer-programming  problem  (BIP). 
We  can  solve  the  BIP  as  a  linear  program  (LP)  (i.e.,  the  linear-programming  relaxation  of 
BIP)  by  ignoring  the  integer  restrictions.  If  the  LP  optimal  solution  contains  only  binary¬ 
valued  variables  we  have  solved  the  BIP.  Otherwise,  we  select  a  fractional  variable,  say 
xk,  and  branch  on  it.  That  is,  we  create  two  new  LPs,  LPi  and  LP2,  by  introducing  the 

constraints  xk  =  0  and  xk  =  1 ,  respectively.  An  optimal  solution  to  the  original  problem  is 

contained  in  the  feasible  region  to  at  least  one  of  these  two  problems.  Next,  we  select  one 
of  the  problems,  say  LPi,  solve  it  (the  LP  relaxation  that  is),  and  evaluate  its  objective.  If 
the  objective  is  no  better  than  the  value  of  some  incumbent  integer  solution  (which  we 
obtained  heuristically,  say),  then  that  problem  is  fathomed  by  bounds.  If  the  problem  is 
infeasible,  it  is  fathomed  by  infeasibility.  If  the  solution  is  integer  then  we  have  fathomed 
by  integrality  (and  we  save  this  solution  as  the  incumbent  if  it  is  better  than  any  we  have 
seen  before).  If  none  of  those  cases  occur  for  LPi,  we  select  a  fractional  variable  and 
recursively  branch  on  it.  Branch  and  bound  stops  when  all  nodes  are  fathomed:  The  best 
solution  found  is  optimal.  (Note  that  some  subproblems  that  have  not  been  solved  may  be 
fathomed  using  information  from  other  subproblem  solutions.  For  instance,  if  an  optimal 
integer  solution  to  LPi  is  found  and  its  objective  value  equals  that  of  LP,  then  LP2  is 
implicitly  fathomed  by  bounding.) 

But  in  our  scheme  for  solving  MXFI,  there  are  no  fractional  values.  Thus,  we  need 
a  different  partitioning  scheme.  We  will  base  the  partition  on,  essentially,  setting  xk  =  0  or 

xk  =  1  for  k  e  At  (A')  ,  where  Ar  (A')  >  R  .  In  this  case,  “forcing  arc  k  into  the  solution”  is 
equivalent  to  setting  xk  =  1  and  “forcing  arc  k  out  of  the  solution”  corresponds  to  setting 
xk  =  0 .  And,  we  will  consider  a  node  fathomed  if,  at  that  node: 

1 .  |  Ar  |  =  R  (equivalent  to  being  fathomed  by  integrality;  an  optimal  solution  to 
the  subproblem  has  been  found),  and 
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2.  z(X,R)  >  UB  .  (Fathomed  by  bounding.) 

If  we  use  this  technique  in  our  problem,  our  enumeration  tree  might  look  like  the 
following,  given  AI(A')  =  {m,  n,  r }  : 


(Level  0) 


(Level  1) 


(Level  2) 


Figure  4.  Sample  Enumeration  Tree  for  R  =  1.  Solutions  are  given  next  to  nodes  in  curly  braces  and 
the  branching  decision  is  next  to  the  branches:  “  The  initial  solution  at  the  root  node  at 
level  0  is  {m,  n,  r} .  “  m  ”  and  “m”  represent  the  exclusion  and  inclusion  of  arc  m, 
respectively.  Note  that  none  of  the  nodes  are  fathomed  in  this  sample  enumeration  tree. 


In  Figure  4,  every  time  we  force  an  arc  in,  we  obtain  the  same  solution.  For 
example,  when  arc  m  is  included  for  the  solution  at  Level  1,  we  get  the  solution  {m,  n,  r } 
from  Level  0.  But,  it  is  unnecessary  to  consider  a  branch  with  {m,  n,  r),  because  we  have 
already  evaluated  this  solution.  So,  we  modify  the  procedure  as  follows: 
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(Level  0) 


(Level  1) 
(Level  2) 


Figure  5.  Modified  Enumeration  Tree  for  R=  1.  Solutions  are  represented  next  to  nodes  in  {}.  mn 
represents  the  inclusion  of  arc  m  and  exclusion  of  arc  n.  A  full  partition  of  the  space  of  all 
possible  solutions  includes  {m,  n,  r},  but  this  branch  need  not  be  created,  because  this 
solution  has  already  been  evaluated  at  the  root  level  of  enumeration  tree. 

In  this  modified  enumeration  tree,  we  begin  with  the  exclusion  of  the  first  arc  m  in 
the  initial  interdiction  set  and  we  obtain  the  solution  {v,  y}  as  a  result.  Since  |  Af  |  A  1  for 
this  solution,  two  new  subproblems  are  introduced  by  adding  the  condition  v  or  vy  .  This 
process  continues  for  all  nodes  until  each  subproblem  node  is  fathomed.  When  all  the 
nodes  are  fathomed,  the  best  feasible  solution  encountered  is  declared  optimal. 

Because  of  limited  time  for  this  study,  we  restrict  the  enumeration  tree  to  a  depth  of 
one  (Level  1).  So,  for  Aj  =  {kvk2,.. .} ,  LB  is  computed  as: 

LB  =  min  {  max  z(A,R\xk  =0),  max  z(A,R  \xk  =  1,  X*  =  0), . . .}  . 

A  *  A  1 

If  the  LB  comes  from  a  node  at  which  \At\  =  R,  At  is  the  optimal  solution  to  MXFI(f?). 
This  is  true  because  all  the  other  nodes  are  fathomed  by  bounding. 

We  find  feasible  solutions  for  each  branch  created  as  follows:  If  |  A{  \  >  R ,  we 
“uninterdict”  arcs  k  e  An  in  order  of  increasing  capacity,  until  \  Ar\  =  R  .  If  |  Aj  \  <  R  ,  we 
interdict  arcs  k  e  {Ac  \  Ar} ,  in  order  of  decreasing  capacity,  until  \  At\  =  R  .  The  feasible 
solution  with  the  lowest  UB  becomes  the  best  feasible  solution. 

If  there  is  more  than  one  missing  value  of  R ,  we  apply  this  procedure  sequentially 
as  follows  ( A+  denotes  the  set  of  arcs  to  be  included): 


28 


Procedure  Branch(  R' ,  R,  A ,  AjR ) ,  UB(r )) 


Input: 


Network  G=  (N,  A),  source  node  5,  sink  node  t, 

R  current  /lvalue  correctly  solved, 

R'  previous  fNvalue  correctly  solved  (R'  >  R  +  1), 

A  value  of  parameter  such  as  A  =  /t*+1  -  a  , 


A\r  1  interdiction  set  that  corresponds  to  x(R') , 

UB(R )  feasible  objective  value  for  x(f?) . 


Output:  x(R) 

A\R) 

LB(R ) 
UB(R) 

absGap(R) 

relGap(R) 


feasible  (possibly  optimal)  solution  for  MXFI(f?), 
interdiction  set  that  corresponds  to  x(R) , 

possibly  improved  lower  bound  on  optimal  objective  value  z(R), 
possibly  improved  feasible  objective  value  for  x(R) , 
for  problematic  R:  absolute  optimality  gap  UB(R)  -  LB(R), 
for  problematic  R:  relative  optimality  gap  100 xabsGap(R)/LB(R). 


{ 

for  (r  =  R'  —  1  to  R  +  1)  {  /*  Solve  for  all  missing  values  of  R  */ 

u'k  <—  min  {uk ,  AR,rk }  for  all  k  e  A ;  /*  adjust  arc  capacities  */ 
LB(r)  <—  00 ; 

A+  <—  0;  /*  This  will  be  the  set  of  arcs  forced  into  the  solution  */ 
for  (a  =  1  to  r  +  1)  {  /*  Include  and  exclude  the  arcs  in  A/(r)  */ 

Let  k  be  the  a'h  arc  in  A/(  r  ); 

If  (  k  <£  A+  )  {  /*  if  k  is  not  already  forced  in  */ 

u[  <r-  uk ;  /*  forces  arc  k  out  of  the  cut,  i.e.,  sets  xk  =  0  */ 

} 

(. Ac,f )  <—  maxFlow  ( G ,  5,  t,  u' ); 

Aj  <—  {k  e  Ac  |  uk=  A}\ 

Aj  <-  Ak  u  A+ ; 

xk  <—  1  Vk  e  Ak;  xk  < —  0  \/k  e  A\Aji 
UBlocal  <—  /- 1  A,  I  Ar,  ; 
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/*  Find  a  feasible  solution  */ 

if  (14  l>r){ 

for  (\Aj\  down  to  r)  { 

^,n,„  <-  ars  mm 11  k ; 


UB local  <-  UBlocal  +  u. 


4  AI~{k min)  5 

xk  <—  0 ; 


} 


} 

/*  Find  a  feasible  solution  */ 

if  (1 4  l<0  { 

for  ( r  down  to  |  Ak  | )  { 


ki 

UB 


<—  arg  max  u. ; 

k^{AcIA,} 


local  UBlocal  Ukmm  ’ 


4  <-  4  +  4max  }  i 
<- 1  ; 

} 

} 

LB local  =f-K'(r~a+\)-, 


LB(r )  =  min  {LB(r),LBlocJ 
UB(r)  =  mm{UB(r),UBtocal} 

u[  =  0;  /*  Set  xk  =  1  in  the  following  subproblems  */ 

A+  <-  A+  u{k}; 


absGap  (r)  <—  UB  (r)  -  LB  (r  ) ; 
relGap  (r)  <—  1  OOxabsGap  (r)/LB  (r) ; 

print  (“  For  problematic  R  =  ”,  r,  “the  approximate  soln.  to  MXFI(f?)  is:  ”); 
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print  (  x ,Aj,UB  (r) ,  LB  (r) ,  absGap  (r) ,  relGap  (r) ); 

} 

} 

In  this  procedure,  we  force  an  arc  out  of  a  solution  by  taking  its  capacity  back  to  its 
original  level,  uk  >  XR, .  This  ensures  that  it  will  never  be  seen  as  interdicted  since  AI  =  {  k 
e  Ac  |  u[  =  XR,}.  “Forcing  in”  is  accomplished  by  setting  u'k  =  0 ,  which  means  that  arc  k  is 
destroyed  and  cannot  be  used. 

Although  it  is  not  shown  in  the  pseudo-code  for  simplicity,  if  \  Ar  \  ^  r ,  we  keep 
solving  max-flow  problems  for  X  e  uk  to  find  \A,(X)\  =  r .  If  it  doesn’t  exist,  then  we 
find  Xmm  and  Xmax  such  that  |  Xf(Xmm)  \<r  <  \AI{Xm!a)\.  After  finding  these  values,  we 
compute  LBlocal  the  way  it  is  computed  in  ProblematicRl(  ): 

LBlocal  =max{(/(Amax)-4ax(^-«  +  1)X(/(^ffii„)-^mi„(^-«  +  l))}  • 

If  we  find  LB  =  UB  ,  then  MXFI  is  solved  exactly  for  that  value  of  R.  Otherwise,  we 
may  obtain  an  improved  LB  and/or  UB ,  which  decreases  the  optimality  gap. 
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IV.  COMPUTATIONAL  RESULTS 


In  this  chapter,  we  describe  the  design  of  test  problems,  and  provide  computational 
results  for  our  Lagrangian-based  procedures  to  solve  MXFI. 

A.  TEST  NETWORK  DESIGN 

We  use  grid  networks  to  test  our  procedures.  “NET  a  xb  ”  denotes  such  a  network, 
where  a  is  the  number  of  nodes  on  the  horizontal  axis  and  b  is  the  number  of  nodes  on  the 
vertical  axis.  The  number  of  nodes  in  a  grid  network  is  |  N  |  =  ab  +  2 ,  and  the  number  of 
arcs  is  \A\  =  a(5b  -  4)  -  5b  +  6  .  A  representation  of  a  small  network  is  given  in  Figure  6: 


Figure  6.  NET  3x3.  Arc  capacities  uk  and  resources  r,  are  not  shown  in  this  network.  The  source 
node  and  sink  node  are  numbered  1  and  2,  respectively. 


As  shown  in  Figure  6,  all  horizontal  arcs  he  in  the  direction  of  “west  to  east,”  and 
diagonal  arcs  are  oriented  in  the  southeast  and  northeast  directions.  There  are  also  vertical 
arcs,  north-pointing  and  south-pointing,  between  adjacent  nodes  in  a  column  of  nodes; 
although  these  are  omitted  in  the  first  and  last  columns,  where  they  would  be  superfluous 
for  maximizing  flow.  Arcs  going  out  of  the  source  node  and  going  into  the  sink  node  have 
infinite  capacity  and  cannot  be  interdicted.  (In  fact,  uninterdictable  arcs  have  rk  =  999 , 


33 


which  effectively  makes  them  uninterdictable.)  All  the  other  arcs  are  interdictable.  Finite 
arc  capacities  uk  are  drawn  from  a  discrete  unifonn  distribution  with  range  [1,50]. 

We  use  the  network  generator  written  in  C  by  Derbes  (1997)  and  modify  it  for 
different  networks.  The  network  files  generated  by  this  program  provide  the  following 
information  about  the  structure  of  the  network:  number  of  nodes,  number  of  arcs,  source 
and  sink  nodes.  Also,  each  arc  k  receives  a  capacity  uk ,  cost  of  repair  ck ,  resource  needed 

to  interdict  rk ,  and  “reparability”  pk .  However,  we  use  only  u,  and  rk  in  our  tests. 

B.  RESULTS 

We  use  six  different  grid  networks  for  testing.  The  smallest  one,  NET  5x5,  has  27 
nodes  and  86  arcs.  The  largest  one,  NET  20x20 ,  has  402  nodes  and  1826  arcs.  All  tests 
are  performed  on  a  533  MHz  Pentium  III  personal  computer  with  192  MB  of  RAM  and 
running  the  Microsoft  Windows  XP  Professional  operating  system.  All  programs  are 
written  in  the  Java  1.3  programming  language. 

We  use  three  different  methods  to  test  the  effectiveness  of  our  Lagrangian-based 
procedures.  Each  method  uses  Algorithm  A  with  a  different  procedure  or  procedures. 
Procedure  ProblematicRl()  is  the  heuristic  that  Bingol  (2001)  uses  to  find  a  feasible 
solution,  and  a  lower  bound  on  z(Xk,R) .  We  use  ProblematicR2()  to  find  the  best 
Lagrangian  multipliers.  The  procedure  Branch()  applies  the  restricted  branch-and-bound 
technique  introduced  in  Chapter  III.  Method  1  uses  ProblematicRl()  integrated  into 
Algorithm  A;  this  is  the  technique  Bingol  (2001)  analyzes.  In  particular,  every  time 
Algorithm  A  runs  into  a  contiguous  set  of  problematic  R -values,  ProblematicRl()  is  called 
to  find  UB  and  LB  for  each.  Method  2  works  the  same  way,  but  ProblematicR2()  is 
invoked  to  find  /t*  for  problematic  R -values  before  parameters  are  passed  on  to 
ProblematicRl()  which  is,  in  this  case,  only  used  to  find  feasible  solution.  Furthermore, 
the  indicated  call  to  Branch()  inside  of  ProblematicR2()  is  not  made.  Method  3  invokes 
procedure  ProblematicR2()  exactly  as  stated  in  the  text,  i.e.,  with  a  call  to  Branch()  if  a 
positive  gap  remains. 
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Since  we  are  really  only  concerned  with  problematic  i? -values,  the  tables  only 
present:  the  problematic  values  R ,  total  run-times,  upper  bounds  z(R),  lower  bounds 
z(A,R),  absolute  gaps  “absGap,”  and  relative  gaps  “relGap.”  Tables  also  include  the 
uninterdicted  maximum  flow  value  /* . 

Table  1  shows  results  for  Method  1.  Although  all  absolute  gaps  are  small,  relative 
gaps  can  be  as  large  as  25.0%  when  the  divisor  z(A,R)  is  small.  Each  network  that  we 
have  tested  has  yielded  one  or  more  problematic  values.  Run-times  are  small,  the  largest 
one  being  31.5  cpu  seconds. 

Table  2  exhibits  results  for  Method  2  using  the  same  networks.  Since  the  best 
Lagrangian  mulitiplier  is  not  necessarily  an  integer,  z(A*r,R )  might  have  fractional  values. 
“absGap”  and  “relGap”  are  the  results  computed  with  these  fractional  values.  However, 
optimal  objective  values  must  be  integral,  so  it  makes  sense  to  round  z(A*r,R )  up  to  the 

nearest  integer.  Thus,  “[absGap]”  and  “[relGap]”  denote  absolute  and  relative  gaps 

computed  after  this  rounding.  Of  36  problematic  R -values  found  with  Method  1,  Method  2 
solves  16  exactly.  The  largest  remaining  relative  gap  is  14.3%  for  R  =  5 1  in  NET  20  x  20 . 
Run-times  increase  only  slightly  for  the  networks  that  have  a  few  problematic  R -values.  If 
ProblematicR2()  is  invoked  many  times,  the  increase  is  about  25%. 

Table  3  presents  the  results  of  Method  3  which  invokes  procedure  Branch().  We 
only  pass  on  to  this  procedure  the  problematic  R -values  remaining  after  Method  2.  With 
this  method,  the  number  of  remaining  problematic  R -values  values  goes  down  to  four  for 
NET  20x20.  The  largest  relative  gap  drops  to  7.1%,  since  a  better  feasible  solution  is 
found  by  Branch().  There  can  be,  unfortunately,  a  drastic  increase  in  run-time;  for 
example,  the  run  time  increases  from  31.5  seconds  to  4013.7  seconds  for  NET  20x20. 

To  explore  the  effectiveness  of  these  methods  further,  we  use  NET  8x15  with 
restricted  capacities.  There  are  four  types  of  restricted  capacities,  denoted  16[  1 ,2],  8[1,4], 
4[1,8],  2[1,16]:  The  notation  “/>[l,g]”  indicates  that  the  random  capacities  are  drawn 
uniformly  from  the  set  {p,2p,...,qp}.  Bingol  (2001)  finds  problem  instances  like  these, 
with  few  distinct  capacities,  to  be  particularly  difficult.  That  is,  they  tend  to  yield  large 
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optimality  gaps.  We  apply  the  same  three  methods  as  before:  Table  4  includes  only 
problematic  /?- values  with  relative  gaps  higher  than  5%  remaining  after  Method  1. 

We  start  with  a  total  of  16  problematic  values  from  Method  1.  The  largest  relative 
gap  from  Method  1  is  20%.  By  using  Method  2,  the  number  of  problematic  /?- values  goes 
down  to  eight,  with  the  largest  gap  being  6.1%.  Four  of  these  problematic  values  are 
solved  exactly  by  Method  3.  However,  there  is  no  improvement  in  the  relative  gaps  of  the 
remaining  four. 
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Table  1.  Test  Results  for  Method  1.  This  is  essentially  Bingol’s  procedure  (Bingol  2001)  which  uses 

Algorithm  A  with  ProblematicRl(),  but  not  ProblematicR2()  or  Branch)).  /*  shows 
maximum  flow  without  interdiction.  The  column  labeled  “R”  gives  the  range  of 
interdiction  resource  levels  tested;  he  maximum  flow  goes  to  0  when  R  equals  the  second  of 
the  two  numbers.  Run  times  cover  the  whole  procedure,  but  only  problematic  /Aval  lies 
remaining  after  the  procedure  is  run  are  shown.  z(R)  is  the  upper  bound  obtained  from  the 
best  solution  with  ProblematicRIQ  and  z(A,R)  is  the  largest  value  of  Lagrangian  function 
obtained  from  MXFI-LRD.  z(A,R)  is  always  integer  since  X  =  w.  for  keA. 
“absOap^’is  z(R)  -  z(A,R)  and  “  relGapj  ”  is  100  x  absGapj  /  z(A,R) . 
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G  =  (TV,  A  ) 

f 

R 

Problematic 

R  -Values 

z(R) 

z(A',R) 

absGap 

relGap 

\z(K,r)\ 

absGap, 

relGap, 

relGap, 

6 

■msa 

m 

1.5 

2.8% 

54 

1 

1.9% 

mmm 

NET  5x  5 

234 

[1,13] 

0.1 

wmttM 

KSUH 

2 

■IS 

■mi 

6.5 

4.3% 

153 

6 

3.9% 

■*aa» 

17 

18 

17.2 

0.8 

4.7% 

18 

- 

- 

5.9% 

16 

29 

27.4 

1.6 

5.8% 

28 

1 

3.6% 

7.4% 

15 

39 

37.6 

1.4 

3.7% 

38 

1 

2.6% 

5.4% 

NET  8x  8 

416 

[1,22] 

0.7 

14 

48 

47.8 

0.2 

48 

- 

- 

2.1% 

6 

202 

200 

2.0 

200 

2 

1.0% 

5 

228 

227 

1.0 

227 

1 

0.4% 

NET  8x12 

15581 

[1,341 

1.0 

15 

189 

188.5 

0.5 

0.3% 

189 

- 

- 

0.5% 

40 

5 

4.3 

0.7 

16.3% 

5 

- 

25.0% 

39 

7 

6.7 

0.3 

4.5% 

7 

- 

16.7% 

34 

27 

26.5 

0.5 

1.9% 

27 

- 

3.8% 

24 

146 

145.5 

0.5 

mmm 

146 

_ 

0.7% 

NET  8x15 

841 

[1,43] 

2.3 

23 

165 

163 

2.0 

1.2% 

163 

2 

1.2% 

1.9% 

22 

181 

180.5 

0.5 

181 

- 

0.6% 

17 

279 

278.5 

0.5 

279 

- 

0.4% 

6 

590 

588.5 

1.5 

589 

0.2% 

0.3% 

38 

11 

10.5 

j&SaKkJHM 

11 

- 

10 

431 

430.3 

nifl 

431 

_ 

NET  15x  15 

878 

[1,43] 

10.3 

9 

463 

462.7 

0.1% 

463 

- 

0.2% 

3 

686 

684.5 

1.5 

0.2% 

685 

0.1% 

0.3% 

53 

7 

7 

- 

- 

7 

- 

- 

14.3% 

51 

16 

14 

2.0 

14.3% 

14 

2 

14.3% 

14.3% 

50 

19 

18 

1.0 

5.6% 

18 

1 

5.6% 

5.6% 

41 

83 

80.5 

2.5 

3.1% 

81 

2 

2.5% 

3.8% 

40 

93 

90 

3.0 

3.3% 

90 

3 

3.3% 

4.5% 

39 

101 

99.5 

1.5 

1.5% 

100 

1 

1.0% 

3.1% 

38 

110 

109 

1.0 

0.9% 

109 

1 

0.9% 

1.9% 

NET  20X20 

1173 

[1,58] 

39.3 

37 

119 

118.5 

0.5 

0.4% 

119 

- 

0.8% 

32 

175 

174.5 

0.5 

0.3% 

175 

- 

0.6% 

17 

508 

507.8 

0.2 

508 

- 

0.2% 

16 

542 

541.5 

0.5 

0.1% 

542 

- 

0.2% 

15 

576 

575.3 

0.7 

0.1% 

576 

- 

0.2% 

14 

609 

609 

- 

- 

609 

- 

0.3% 

11 

712 

711.5 

0.5 

0.1% 

712 

- 

0.1% 

7 

859 

858 

1.0 

0.1% 

858 

0.1% 

0.1% 

Table  2. 


Test  Results  for  Method  2.  This  method  uses  Algorithm  A  and  invokes  ProblematicR2() 

* 

and  ProblematicRIQ.  ProblematicR2()  finds  the  best  Lagrangian  multipliers  AR  for 
problematic  R -values,  which  are  then  passed  on  to  ProblematicRIQ  to  find  feasible 


solutions.  Since  z(R)  is  integer, 


is  a  valid  lower  bound  on  z(R).  “  absGap2  ”  is 


z(R)~ 


<\,R) 


and  “re!Gap2”  is  100  x  absGap,/ 


z(\,R) 


‘  absGap ' 


is 


z(R)-z(/ 1  ,R)  and  “relGap”  is  100 x absGap / z(A,  ,R).  indicates  that  MXFI  is 

R  R 

solved  optimally.  Only  two  problematic  /Aval lies,  R  =  53  and  R  =  14  for  NET  20  x  20  ,  are 

* 

solved  exactly  before  rounding  up  z(An,R)  but,  overall,  more  than  50%  of  the  problematic 

R 

A'- values  are  solved  with  rounding.  relGap  j  is  shown  for  comparison. 
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G  =  (N,  A) 

* 

R 

Run  Time 

Problematic 

\LB-\ 

absGapj 

relGap  3 

relGap, 

relGap, 

/ 

(cpu  sec.) 

R  -Values 

z(R) 

LB 

NET  5X5 

234 

[1,13] 

0.2 

6 

55 

55 

55 

- 

- 

3.8% 

1.9% 

2 

159 

158 

158 

1 

0.6% 

4.6% 

3.9% 

17 

18 

18 

18 

- 

- 

5.9% 

- 

16 

29 

29 

29 

- 

- 

7.4% 

3.6% 

15 

39 

39 

39 

- 

- 

5.4% 

2.6% 

NET  8x8 

416 

[1,22] 

7.0 

14 

48 

48 

48 

- 

- 

2.1% 

- 

6 

202 

202 

202 

- 

- 

1.0% 

1.0% 

5 

228 

228 

228 

- 

- 

0.4% 

0.4% 

NET  8x12 

690 

[1,341 

4.4 

15 

189 

189 

189 

- 

- 

0.5% 

- 

40 

5 

4.3 

5 

- 

- 

25.0% 

- 

39 

7 

6.7 

7 

- 

- 

16.7% 

- 

34 

27 

27 

27 

- 

- 

3.8% 

- 

NET  8x15 

841 

[1,43] 

25.2 

24 

146 

146 

146 

- 

- 

0.7% 

- 

23 

165 

165 

165 

- 

- 

1.9% 

1.2% 

22 

181 

181 

181 

- 

- 

0.6% 

- 

17 

279 

279 

279 

- 

- 

0.4% 

- 

6 

590 

590 

590 

- 

- 

0.3% 

0.2% 

38 

11 

11 

11 

- 

- 

mm 

- 

NET  15x15 

878 

[1,43] 

96.0 

10 

431 

431 

431 

- 

- 

Kl 

- 

9 

463 

463 

463 

- 

- 

■n 

- 

3 

686 

686 

686 

- 

- 

0.3% 

0.1% 

51 

15 

14 

14 

1 

14.3% 

14.3% 

50 

19 

18 

18 

1 

5.6% 

5.6% 

5.6% 

41 

83 

83 

83 

- 

- 

3.8% 

2.5% 

40 

93 

92 

92 

1 

1.1% 

4.5% 

3.3% 

39 

101 

101 

101 

- 

- 

3.1% 

1.0% 

38 

110 

110 

110 

- 

- 

1.9% 

0.9% 

NET  20x20 

1173 

[1,58] 

4013.7 

37 

119 

119 

119 

- 

- 

0.8% 

- 

32 

175 

174.5 

175 

- 

- 

0.6% 

- 

17 

508 

508 

508 

- 

- 

0.2% 

- 

16 

542 

542 

542 

- 

- 

0.2% 

- 

15 

576 

576 

576 

- 

- 

0.2% 

- 

11 

712 

712 

712 

- 

- 

0.1% 

- 

7 

859 

859 

859 

- 

- 

0.1% 

0.1% 

Table  3.  Test  Results  from  Method  3.  Method  3  is  essentially  Method  2  except  that  Branch()  is 
invoked  for  problematic  R-values  that  Method  2  cannot  resolve.  LB  =  min 
{max  z(A, R  \  xk  =  0),  max  z(A,R  \  xk  =  1 ,  =0),...}.  “absGap3”is  z(2)  -  \ LB~\  and 

“relGap3”  is  100xabsGap3  /  [ LB~\ .  Absolute  and  relative  gaps  computed  with  fractional 

values  are  not  shown  in  this  table.  12  problematic  R -values  out  of  the  16  remaining  after 
Method  2  are  solved  exactly  by  Method  3.  UB  for  R  =  51  of  NET  20  x  20  drops  from  16 

to  15,  because  a  better  feasible  solution  is  found  by  BranchQ.  reKiap  |  and  relGap2  are 
also  shown  in  the  table  for  comparison. 
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Method  1 

Method  2 

Method  3 

G=(N,  A ) 

Problematic 

R  -Values 

z(R) 

z(A,R) 

relGapj 

z(R) 

\<A"r,R)\ 

relGap, 

z(R) 

r^i 

relGap, 

24 

320 

304 

5.3% 

304 

304 

- 

- 

- 

- 

23 

352 

320 

10.0% 

320 

320 

- 

- 

- 

- 

22 

384 

336 

14.3% 

336 

336 

- 

- 

- 

- 

1 6[  1 ,2] 

21 

400 

352 

13.6% 

352 

352 

- 

- 

- 

- 

20 

416 

368 

13.0% 

368 

368 

- 

- 

- 

- 

19 

432 

384 

12.5% 

400 

394 

1.5% 

400 

394 

1.5% 

18 

448 

400 

12.0% 

432 

420 

2.9% 

432 

420 

2.9% 

17 

464 

432 

7.4% 

464 

446 

4.0% 

464 

446 

4.0% 

31 

104 

96 

8.3% 

96 

96 

- 

- 

- 

- 

8[1,4] 

30 

120 

104 

15.4% 

112 

110 

1.8% 

112 

112 

- 

29 

128 

120 

6.7% 

128 

123 

4.1% 

128 

123 

4.1% 

4[ 1,8] 

38 

240 

200 

20.0% 

200 

200 

- 

- 

- 

- 

37 

280 

240 

16.7% 

240 

240 

- 

- 

- 

- 

35 

280 

260 

7.7% 

280 

264 

6.1% 

280 

280 

- 

2[1,16] 

34 

340 

320 

6.3% 

340 

329 

3.3% 

340 

340 

- 

33 

400 

380 

5.3% 

400 

393 

1.8% 

400 

400 

- 

Table  4.  Comparative  Test  Results  for  NET  8x15  with  Restricted  Arc  Capacities.  “Methods”  are  defined  in  the  text.  Arc  capacities  denoted  p[l,q] 

are  drawn  from  the  set  {  p,2p,...,cjp  }.  “relGapj  ,relGap2,  and  relGap3  ”  are  the  relative  gaps  obtained  from  Method  1,  2  and  3,  respectively. 

Only  R-values  with  relative  gaps  remaining  greater  than  5%  after  Method  1  are  shown  in  this  table.  50%  of  problematic  R-values  are  solved 
exactly  by  Method  2,  and  50%  of  those  remaining  are  then  solved  by  Method  3.  Notice  that  UB  for  R  =  30  of  the  network  with  arc  capacities 
8[1,4]  has  dropped  from  120  to  112,  because  R=3l  is  solved  exactly  by  Method  2,  and  this  makes  it  possible  to  find  a  better  feasible  solution 
for  R  =  30. 
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V.  CONCLUSIONS  AND  FURTHER  RESEARCH 


A.  SUMMARY 

The  maximum-flow  network-interdiction  problem  (MXFI)  arises  when  an 
“interdictor”  uses  limited  interdiction  resources  to  restrict  an  adversary’s  use  of  a  network. 
Conceptually,  the  interdictor  destroys  arcs  to  minimize  the  maximum  flow  that  the 
adversary,  or  “network  user,”  can  achieve  through  the  network.  In  this  study,  we  have 
explored  new  Lagrangian-relaxation  techniques  for  solving  MXFI  for  all  values  of 
interdiction  resource  R  in  a  specified  range.  We  assume  that  one  unit  of  resource  is 
required  to  interdict  an  arc,  so  we  only  deal  with  integer  values  of  R. 

MXFI  can  be  solved  by  means  of  integer  programming,  but  this  approach  is  slow. 
Bingol  (2001)  modifies  the  Lagrangian-relaxation  technique  developed  by  Derbes  (1997) 
to  solve  MXFI  for  all  R  e  [1 , |  Amm  |],  where  Amm  denotes  a  minimum  cardinality  cut 

consisting  of  only  interdictable  arcs.  The  Lagrangian-relaxation  technique  solves  a 
sequence  of  max-flow  problems  with  related  arc  capacities,  and  is  therefore  quite  efficient. 

Let  z(R )  denote  the  optimal  objective  value  to  MXFI  for  a  given  R,  and  let 
z{A,R)  denote  the  Lagrangian  lower  bound  on  z(R)  given  Lagrangian  multiplier  A. 
Essentially,  Bingol  (2001)  evaluates  z(A,R)  simultaneously  for  all  R,  but  only  for  values 
of  A  equal  to  perturbed  arc  capacities.  We  follow  this  procedure — the  range  of  A  is 
correct — but,  when  z{A,R)  <~z(R)  for  A*=  argmax  z(A,R) ,  where  z(R )  is  the  objective- 

A 

function  value  for  the  best  feasible  solution  found.  This  is  accomplished  with  a  search 
technique  based  on  the  slope  of  z{A,R)  with  respect  to  A  .  For  most  test  problems,  we  can 
decrease  the  optimality  gap  for  “problematic  /Lvalues”  by  approximately  50%.  A  value  of 
R  is  “problematic”  if  Bingol’s  procedure  leaves  z(A,R)  <z(R).  Although  this  approach 
focuses  on  tightening  lower  bounds,  it  is  also  finds  better  feasible  solutions  on  occasion, 
and  sometimes  proves  optimality  of  the  best  solution.  However,  some  large  optimality 
gaps  may  remain. 
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We  tackle  the  remaining  problematic  values  with  a  modified  branch- and-bound 
procedure.  This  procedure  uses  a  partition  based  on  the  best  interdiction  set  found, 
Aj  ={kvk2,...,kL) .  In  particular,  we  solve  max  z(A,R  \  xk  =  \,xk  =  0),..., 

A  1 

max  z(A,  R  \  xki  =  \,xk^  =1,  . . . ,  xk/ 1  =  1,  xk/  =  0) .  ( max  z(A,  R  \  xk  =  1  \/k  e  Ak )  is  already 

solved.)  The  minimum  lower  bound  from  these  problems  is  a  lower  bound  on  the  optimal 
solution.  This  partitioning  scheme  can  be  applied  recursively  to  solve  MXFI  exactly  for  a 
given  R,  but  we  have  only  implemented  a  single  application  of  the  partition.  That  is,  we 
have  limited  the  enumeration  tree  of  this  method  to  one  level  below  the  root.  Using  this 
technique,  we  may  find  a  lower  bound  which  equals  the  upper  bound  obtained  from  the 
best  feasible  solution.  In  this  case,  MXFI  is  solved.  In  our  tests,  we  have  solved 
approximately  70%  of  the  problematic  /U values  that  remain  after  applying  the  first 
technique.  But,  this  procedure  runs  slowly  because  of  the  need  to  solve  many  max-flow 
problems,  so  it  is  best  used  only  when  large  optimality  gaps  are  encountered. 

Having  used  these  two  techniques  together,  we  have  been  able  to  solve  MXFI  either 
exactly  or  with  small  optimality  gaps  for  most  of  the  problematic  R -values  encountered. 
But,  we  still  cannot  guarantee  an  exact  solution  or  a  small  optimality  gap. 

B.  FURTHER  RESEARCH 

In  this  thesis,  we  have  used  a  restricted  modified  branch-and-bound  procedure  that 
explores  the  enumeration  tree  only  one  level  below  the  root.  That  is  why  some  problematic 
f? -values  still  remain.  A  fully  recursive  implementation  of  the  branch-and-bound  algorithm 
will  solve  MXFI  exactly,  at  least  in  theory.  Even  the  restricted  branch-and-bound 
procedure  becomes  slow  because  of  the  need  to  solve  many  max-flow  problems.  Thus,  a 
fully  recursive  procedure  based  on  this  methodology  would  be  even  slower.  To  make  such 
a  procedure  reasonably  efficient,  a  faster  maximum-flow  algorithm  should  be  used  and 
techniques  for  reoptimizing  related  maximum-flow  problems  should  be  explored. 

For  the  sake  of  simplicity,  we  have  assumed  that  that  rk  =  1  for  all  interdictable  arcs 
k.  That  is,  one  unit  of  resource  is  needed  to  interdict  any  interdictable  arc.  Allowing 
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rk  >  1  (but  still  integer)  will  make  the  problem  harder.  In  this  case,  a  complete  branch-and- 
bound  procedure  would  become  more  important  to  possess. 

We  have  assumed  that  grid  networks  are  adequate  representatives  for  real-world 
networks,  so  we  have  performed  our  tests  on  these  networks.  However,  we  are  uncertain 
how  well  these  procedures  will  perform  on  actual  networks.  Obviously,  future  research 
should  apply  these  techniques  to  more  realistic  networks. 
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